clc,clear,close all
load('pic1_3_cz.mat');
% % NAME='Pic1_3.bmp';
% % z=sub_pixel(NAME);
i=imadjust(z,[0.2;0.9],[0.9;1],25);%亮度变换
%原图中没有噪点，问题在于高光，所以不需要滤波
i=imbinarize(i);
i1=imclearborder(i);%提取高光部分
i=i-i1;
i=imadjust(i,[0;1],[1;0]);
% imshow(i)
% se = strel('square',3);
% se = strel([0 1 0;1 1 1;0 1 0]);
se = strel('square',30);
for n=1:5
i= imdilate(i,se);
end
 for n=1:5
i= imerode(i,se);
 end
i=imadjust(i,[0;1],[1;0]);
a=edge(i,'sobel');%使用sobel算子进行边缘检测
figure;
imshow(a);
